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Low  frequency  Av  =  0-350  cm"1,  Raman  intensity  data  were  obtained  from  liquid  water 
between  3.5  and  89  3  °C  using  holographic  grating  double  and  triple  monochromators.  The 
spectra  were  Bose-Einstein  (BE)  corrected,  7/(1  +  n),  and  the  total  integrated  (absolute) 
contour  intensities  were  treated  by  an  elaboration  of  the  Young-Westerdahl  (YW) 
thermodynamic  model,  assuming  conservation  of  hydrogen-bonded  (HB)  and 
nonhydrogen-bonded  (NHB  =  bent  and/or  stretched,  O-H  O)  nearest-neighbor  0-0  pairs.  A 
A H\  value  of  2.6  +  0.1  kcal/mol  O-H  —  O  or  5.2  +  0.2  kcal/mol  H:0  ( 1 1  kJ/mol  O-H  —  O, 
or  22  kJ/mol  H20)  resulted  for  the  HB— »NHB  process.  This  intermolecular  value  agrees 
quantitatively  with  Raman  and  infrared  A H°  values  from  one-  and  two-phonon  OH-stretching 
regions,  and  from  molecular  dynamics,  depolarized  light  scattering,  neutron  scattering,  and 
ultrasonic  absorption,  thus  indicating  a  common  process.  A  population  involving  partial 
covalency  of,  i.e.,  charge  transfer  into,  the  H  O  units  of  linear  and/or  weakly  bent  hydrogen 
bonds,  O-H  •••  O;  is  transformed  into  a  second  high  energy  population  involving  bent,  e  g.,  150° 
or  less,  and/or  stretched,  e.g.,  3.2  A,  but  otherwise  strongly  cohesive  O-H  O  interactions-  All 
difference  spectra  from  the  fundamental  OH-stretching  contours  cross  at  th  tX(ZJ(  +  Z)Y 
isobestic  frequency  of  3425  cm'1.  Also,  total  integrated  Raman  intensity  decreases  occurring 
below  3425  cm  - 1  with  temperature  rise  were  found  to  be  proportional  to  the  total  integrated 
intensity  increases  above  3425  cm- ',  indicating  conservation  among  the  HB  and  NHB  OH- 
stretching  classes.  From  the  enthalpy  of  vaporization  of  water  at  0  °C,  and  the  A H  \  of  2.6 
kcal/mol  O-H  •••  O,  the  additional  enthalpy,  A/f  2,  needed  for  the  complete  separation  of  the 
NHB  0-0  nearest  neighbors  is  —  3.2  kcal/mol  O^-H  —  O  or  ~  6.4  kcal/mol  H.O  ( 1 3  kJ/mol 
O-H  —  O  or  26  kJ/mol  H.O).  The  NHB  0-0  nearest  neighbors  are  held  by  forces  other  than 
those  involving  H  O  partial  covalency,  i.e.,  electrostatic  (multipole),  induction,  and 
dispersion  forces.  The  NHB  0-0  pairs  do  not  appear  to  produce  significant  intermolecular 
Raman  intensity  because  they  lack  H  — O  bond  polarizability,  but  the  corresponding  NHB  OH 
oscillators  do  contribute  weakened  Raman  intensity  above  3425  cm  - '.  An  ideal  solution 
thermodynamic  treatment  employing  A H\  —  2.6  kcal/mol  O-H  —  O,  the  HB  mole  fraction, 
and  the  vapor  heat  capacity,  yielded  a  very  satisfactory  specific  heat  value  of  1 . 1  cal  deg  - 1  g~ 1 
H;o  at  0  °C.  The  NHB  mole  fraction, from  the  YW  treatment  is  negligibly  small,  0.06  or 
less,  for  r<  —  50  °C.  However,/,  increases  to  0.16  at  O’C,  and/,  1  at  1437  °C,  where  recent 
shock-wave  Raman  measurements  indicate  loss  of  all  partially  covalent,  charge  transfer 
hydrogen  bonding. 


I.  INTRODUCTION 

The  low-frequency  Raman  spectrum  from  liquid  water, 
Av  —  0-350  cm'1,  is  composed  of  (at  least)  two,  weak, 
broad  features  at  nominal  positions  of  170  and  60  cm'  '.  The 
170  cm" '  band  was  discovered  by  Segre  in  1931'  and  the 
weaker  60  cm" 1  band  by  Bolla  in  1932. 2  Microdensitometer 
tracings  presented  by  Segre  clearly  indicate  that  both  the  1 70 
cm  1  intensity,  and  the  Raman  intensity  near  3250  cm  ~ 1  in 
the  OH-stretching  region,  decrease  relative  to  other  Raman 
bands  with  temperature  rise.  Such  results  were  interpreted 
by  Segre  as  a  breakdown  of  intermolecular  structure,  and 
agree  with  data  of  this  work. 
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Similar  temperature  effects  were  observed  by  Bolla  for 
the  60  cm'  1  band,  as  well  as  for  the  Raman  intermolecular 
librational  and  OH-stretching  bands.2  Because  of  the  high 
quality  of  Bolla’s  low-frequency  spectra,  correction  for  the 
Bose-Einstein  (BE)  thermal  population  factor  and  stray- 
light  background  was  accomplished  recently. '  These  correc¬ 
tions  yielded  agreement  with  modern  results. 

The  1 70  cm  " 1  Raman  band  arises  from  the  vibration  of 
an  oxygen  atom  against  its  nearest-neighbor  oxygen  when  a 
linear  or  at  most  weakly  bent  hydrogen  bond,  O-H  —  O,  is 
present  between  them.  This  motion  is  often  described  alter¬ 
natively  as  a  restricted  translation,  but  one  that  is  perpendic¬ 
ular  to  the  hydrogen-bond  O-H  O  direction,  i.e.,  as  a 
bending  of  the  hydrogen  bond. 

Raman  activity  at  170  and  60  cm" 1  occurs  because  of 
the  partial  covalency  of,4  or  the  charge  transfer  into.'  the 
H  —  O  part  of  the  O-H  —  O  unit.  Without  this  partially  co- 


valent  character,  no  change  of  the  H  -  O  bond  polarizability 
could  occur  during  vibration. 

Quantum  mechanical  calculations,  however,  indicate 
that  the  partially  covalent  character  of  the  hydrogen  bond 
decreases  with  decrease  of  the  O-H  —  O  angle  (below 
180°).'’  This  suggests7  that  the  Raman  polarizability  deriva¬ 
tive  is  a  function  of.  and  decreases  very  sharply  for,  O- 
H  •  O  angles  below  some  cut-off  value,  e.g.,  roughly  150°, 
see  subsequent  discussion. 

In  this  work  linear  or  weakly  bent  hydrogen  bonds  are 
designated  by  O-H  —  O,  or  by  HB  (hydrogen  bonded). 
Strongly  bent  and/or  stretched  (ruptured)  interactions  are 
designated  as  nonhydrogen-bonded.  NHB,  or  by  O-H  O, 
without  dots.  (The  absence  of  dots  indicates  lack  of  partial 
covalency  or  charge  transfer). 

Numerous  authors  have  now  reported  Raman  and  in¬ 
frared  results  for  the  low-frequency  region  of  water. IU2ft 
Moreover,  a  recent  trend  has  been  to  report  frequency  ic- 
dueed  spectra^/ /( 1  n),  where  (1  +  n  )  is  the  one-phonon 

B'  e-Einstein  (BE)  factor  with  n  =  l/(e*“"*r  —  D.rjisthe 
frequency,  and  /  the  Raman  amplitude  at  tu2’-29;  or,  tem¬ 
perature  reduced  spectra/ /( 1  ■+■  n),  where  only  the  BE  cor¬ 
rection  has  been  made. '  These  reduced  spectra  involve  dif¬ 
ferent  functionalities  of  the  frequency  dependent  Raman 
coupling  coefficient  r(ai), 27-29  but  in  many  cases  no  clear 
basis  for  the  assumed  functionality  exists,  because  neutron 
inelastic  scattering  spectra  are  required  to  determine  c(co), 
and  such  spectra  are  often  unavailable.  However,  a  safe  ap¬ 
proach.  employed  here,  is  to  ignore  the  Shuker-Gammon 
formalism  entirely.”  and  simply  to  apply  the  BE  correction, 
as  clearly  indicated  by  Placzek  in  1934.'°  In  any  case,  it 
should  be  made  clear  that  either  the  «//( 1  -+-  n)  or  the 
/ <  1  -  n )  representation  is.  in  principle,  equivalent,  with  re¬ 
gard  to  any  additional  temperature  dependence  that  remains 
after  correction.  The  to  term  is  simply  a  temperature-inde¬ 
pendent  multiplier,  and  therefore  it  cannot  affect  the  cor¬ 
rected  temperature  dependence.  On  the  other  hand,  failure 
to  apply  any  temperature  correction  whatsoever  causes  the 
low-frequency  Raman  intensities  from  water  to  appear  to 
decrease  too  rapidly  with  temperature  rise.12 

In  previous  ( 1966)  low-frequency  work  involving  wa¬ 
ter.  a  XH  °  value  of  about  5.6  keal/mol  was  obtained  from  the 
temperature  dependence  of  the  Raman  intensity.12  BE  cor¬ 
rection  was  attempted,  but  it  was  inadequate,  and  the  result¬ 
ing  XH '  value  was  too  large,  as  seen  from  the  calculated  heat 
capacity.  However,  with  modem  methods  of  data  acquisi¬ 
tion  and  computerized  scanning  and  processing,  it  is  easy  to 
obtain  BE  corrected  spectra,  thus  greatly  reducing  the  labor 
involved  when  many  corrected  spectra  at  a  series  of  tempera¬ 
tures  are  required.  Moreover,  BE  reduction  of  the  total  scat¬ 
tering  amplitude  leads  to  a  surprisingly  simple  and  accurate 
base  line  not  recognized  previously12  in  which  a  strongly 
curving  uncorrected  base  line  was  used.  ( Furthermore,  no 
simple  method  of  obtaining  the  base  line  for  frequency  re¬ 
duced  spectra  was  found  here. )  Thus,  although  the  uncor¬ 
rected  Raman  spectra  obtained  now  and  previously 12  do  not 
differ  substantially,  the  present  BE  correction  yields  accu¬ 
rate  results  which  constitute  the  principal  subject  of  this  new 
work. 


Another  aspect  of  the  present  work  involves  difference 
Raman  spectra  obtained  from  the  fundamental  (one- 
phonon)  intramolecular  OH-stretching  region,  as  well  as 
the  isosbestic  frequency  from  the  X(ZJC  +  Z)Y  spectra. 
Precise  isosbestic  points  were  also  measured  for  the 
X(ZZ)Y,X{  YJ(  +  Z)  YtindX(ZX)  Y orientations.  In  addi¬ 
tion,  precise  isosbestic  points  were  observed  for  /( VV )  —  4/ 
3/(VH),  '  t.,  for  a2,  but  only  the  X(ZJl  •+ -  Z)Y  data  are 
presented  here.  The  difference  method  yields  results  which 
are  entirely  consistent  with  previous  Gaussian  computer  de¬ 
compositions  of  the  OH-stretching  contours.1213 17 19  31  but 
it  does  not  require  assumptions  related  to  component  shapes 
or  to  the  number  of  components. 

The  qualitative  OH-stretchmg,  OH-difference,  and  low- 
frequency  Raman  data  are  now  presented. 

II.  EXPERIMENTAL  PROCEDURES 

Low-frequency  Raman  spectra  were  obtained  with  a 
J.  Y.  T800  triple  monochromator,  and  with  a  J.  Y.  HG2S  ho¬ 
lographic  grating  double  monochromator.  The  T800  is  a  ho¬ 
lographic  grating  instrument  of  800  mm  focal  length.  It  em¬ 
ployes  plane  gratings  in  each  of  the  three  Czemy-Tumer 
type  monochromators  in  series.  With  the  T800,  data  collec¬ 
tion  and  processing  was  accomplished  using  a  pulse  discri¬ 
minator,  a  custom  interface,  a  Commodore  64  computer 
plus  1541  disk  drive,  and  a  HP  7470  plotter.  The  J.Y.  HG2S 
system  employed  an  EG&G  model  1182  preamplifier  discri¬ 
minator,  an  ISA  Spectra  link  photon  counter,  an  ISA  Co¬ 
lumbia  Data  Products  computer,  and  the  HP  7470  plotter, 
although  dc  detection  was  also  accomplished  with  a  Keithly 
4 14S  picoammeter.  BE  corrected  spectra  /  /  ( 1  +  n )  and  fre¬ 
quency  reduced  spectra  tol /( 1  +  n)  were  obtained  with  the 
computers  described. 

Raman  spectra  were  obtained  in  the  low-frequency  and 
OH-stretching  regions  under  rigorously  constant  conditions 
of  excitation  intensity  and  detection  geometry.  These  spectra 
were  thus  absolute,  except  that  they  were  not  calibrated  in 
terms  of  cm2  mol -  1  sr  “ 1  vs  temperature.  This  intensity  con¬ 
stancy  was  verified  by  a  series  of  scanning  cycles  in  which  the 
peak  amplitude,  e  g.,  of  the  OH-stretching  peak,  was  deter¬ 
mined  to  be  the  same  within  39c,  or  less,  when  the  water  was 
returned  to  the  same  temperature  after  a  long  time  interval. 
This  procedure,  which  monitors  long  term  drift,  was  used 
previously,  where  Raman  OH-stretching  overtone  intensi¬ 
ties  were  held  constant  to  within  3%  over  periods  of  10  h.31 

A  special  Raman  cell  was  used  for  measurement  of  the 
isosbestic  frequency  in  the  OH-stretching  region.  This  cell 
was  surrounded  by  an  outer  container  having  two  open  la¬ 
ser-beam  ports  (3  mm  diameter)  and  one  large,  10 X  10  cm, 
glass  viewing  window.  The  outer  container  was  flushed  with 
dry  Ar  to  prevent  condensation  of  moisture  on  the  window  s 
of  the  Raman  cell.  The  Raman  cell  employed  three  windows, 
two  in-line  windows  for  laser  excitation,  and  a  third  large  ( 2 
cm  diameter)  viewing  window  exactly  at  90*  to  the  other 
two.  The  in-line  windows  were  exactly  parallel  to  each  other 
because  they  rested  directly  on  precisely  machined,  parallel, 
stainless  steel  faces  of  the  Raman  cell. 

The  laser-beam  position  in  the  filled  Raman  cel!  was 


held  precisily  fixed  for  the  duration  of  a  series  of  measure¬ 
ments.  Invariance  of  the  laser-beam  position  in  the  water  is 
extremely  important  for  quantitative  intensity  measure¬ 
ments  under  conditions  of  changing  temperature.  This  in¬ 
variance  was  accomplished  by  observing  Newton’s  rings  due  I 
to  the  backward  laser  reflection  from  the  entrance  (side) 
window.  ( Newton’s  rings  were  observed  around  a  pinhole  in 
a  black  card  placed  3  m  from  the  side  window.)  No  cell 
adjustment  whatsoever  was  required,  or  made,  during  the 
measurements.  Moreover,  when  Newton’s  rings  were  invar¬ 
iant.  the  laser  spot  on  the  wall  3  m  from  the  exit  window  was 
exactly  fixed.  If  the  laser-beam  position  in  the  filled  Raman 
cell  is  not  exactly  perpendicular  to  the  entrance  (side)  win¬ 
dow  of  the  cell,  bending  of  the  beam  will  occur  as  the  refrac¬ 
tive  index  of  the  water  changes  with  temperature,  and  large 
intensity  errors  will  result.  Such  errors  can  completely  pre¬ 
clude  the  observation  of  precise  isosbestic  points,  and  may 
have  contributed  to  an  erroneous  conclusion  about  the  ab¬ 
sence  of  the  a*  isosbestic  point  for  water.”  The  temperature 
was  maintained  to  ±0.1  ‘C  by  pumping  water  or  cold  brine 
from  an  external  themostated  bath  through  the  large  rear 
chamber  of  the  cell.  However,  the  temperature  of  the  sample 
water  was  measured  accurately  by  placing  the  tip  of  a  ther¬ 
mocouple  1  mm  above  the  laser  beam. 

The  water  employed  was  distilled  three  times  in  a  fused 
quartz  still.  This  triple  distilled  water  was  found  to  be  free  of 
particulates,  and  to  have  an  extremely  low  fluorescence  lev¬ 
el.31 

A  du  Pont  310  analog  computer  was  used  for  Gaussian 
analysis  of  low-frequency  Raman  contours. 

HI.  EXPERIMENTAL  RESULTS,  DISCUSSION,  AND 
INTERPRETATION 

A.  Low-frequency  Raman  data 

I 

Temperature  reduced  or  BE  corrected  Raman  spectra 
from  water,  //( 1  +  n),  obtained  from  —0-500 cm-1  under 
rigorously  constant  conditions  of  excitation  and  detection, 
i.e..  absolute  in  the  sense  described  in  Sec.  II.  are  shown  in 
Fig.  1.  Intensity  maxima  near  53  and  175  cm'1,  i.e.,  the 
nominal  60  and  1 70  cm  " 1  features,  are  obvious  in  the  spec¬ 
trum  corresponding  to  3.5  ’C,  top  spectrum.  A  marked  de¬ 
cline  in  the  intensity  of  the  1 70  cm  " 1  feature  is  also  evident 
with  temperature  rise  to  89.3  "C,  bottom  spectrum. 

Figure  1  shows  four  spectra  in  the  temperature  range 
from  3.5  to  89.3  'C  in  which  a  total  of  12  quantitative  tem¬ 
perature-reduced  spectra  were  obtained.  A  single  constant 
average  base  line  was  determined  from  examination  of  these 
1 2  spectra,  and  tested  as  described  next,  see  the  curved  lower 
lines  in  Fig.  1 .  This  constant  base  line  was  drawn  under  all  of 
the  spectra  by  identical  positioning  of  a  metal  template, 
shaped  specifically  to  reproduce  the  average  baseline  exact¬ 
ly  The  area  above  this  base  line  from  0  to  —325  cm"1 
constitutes  the  total  integrated  contour  intensity,  whose 
temperature  dependence  refers  to  the  HB—  NHB  transfor¬ 
mation  in  water. 

The  base  line  shown  in  Fig.  1  arises  from  the  stray-light 
level  of  the  T800  triple  monochromator.  This  base  line  was 
determined  and  tested  by  comparisons  between  the  low-fre- 
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FIG  J  Bose-Einstein  corrected  low-frequency  Raman  (absolute) 
XiZ.Z  -  X)  >' spectra  from  water  between  3.5  and  8^  3  °C  under  rigorously 
constant  conditions  of  excitation  intensity  and  detection  geometry  Note 
that  because  of  superposition,  the  photon  counts  per  second  are  not 
expressed  at  800  c/s  for  the  top  three  spectra. 


quency  BE  corrected  Raman  spectra  from  water  and  from 
96  wt.  %  H-SOj,  both  at  22  °C.  and  under  the  same  excita¬ 
tion  and  detection  conditions.  The  base  line  was  tested 
further  by  obtaining  low-frequency  BE  corrected  Raman 
spectra  from  water  using  the  HG2S  double  monochromator. 

The  low-frequency  BE  corrected  Raman  spectrum  from 
96%  H;S04  displays.  ( 1 )  a  pronounced  asymmetric  feature 
peaking  sharply  at  25  cm"1  and  showing  a  long  region  of 
downward  concavity  extending  to  150  cm-1,  (2)  a  very 
weak  (compared  to  water),  broad  O-H  —  O  stretching  con¬ 
tour  centered  near  155  ±  5  cm" \  (3)  an  inflection  at  115 
cm'1  between  the 25 and  155 cm" 1  peaks, and  (4)  an  inten¬ 
sity  minimum  at  280  cm" '.  All  of  these  features  occur  below 
the  H:S04  and  HSO/  deformation  frequencies  of  395  and 
419  cm" '  (unresolved).  The  sharply  peaking  nature  of  the 
25  cm”'  feature  allowed  the  total  scattering  amplitude  near 
Av  =  0  cm" 1  to  be  determined.  The  total  scattering  ampli- 


TABLE  I.  Bc»e- Einstein  corrected  low-frequency  total  integrated  contour 
intensities.  /,  for  water  at  various  temperatures.  /*C. 


t.'C 

I 

3.5 

21.46 

11.2 

23.51 

16.7 

23.68 

22.2 

22.17 

300 

22.35 

34.3 

22.23 

464 

17,61 

530 

19.32 

61.8 

18.99 

69.5 

19.91 

78.3 

16.99 

89  3 

18.16 

tude  at  280  cm-  '  could  be  measured  directly.  Further,  the 
base  line  at  the  115  cm"'  inflection  could  be  estimated  be¬ 
cause  of  the  shape  and  weakness  of  the  neighboring  155  +  5 
cm-1  O-H  —  O  feature.  The  entire  base  line  from  96% 
H:S04  was  thus  determined  as  a  smooth  curve  having  the 
same  shape  and  upward  concavity  as  that  shown  in  Fig.  1. 
Furthermore,  the  relative  amplitude  of  the  96%  H,S04  base 
line  was  found  to  match  that  obtained  for  water  at  22  °C  to 
within  +  5%  over  the  entire  range  0-300  cm  ~ 

The  base  line  slope  of  Fig.  1  was  also  tested  with  the 
HG2S  double  monochromator.  For  water,  the  same  shape 
and  upward  concavity  was  observed  w  ith  both  of  the  HG2S 
and  T800  instruments.  The  ratio  of  the  total  scattering  am¬ 
plitude  near  Av  =  0  to  that  at  A  v  =  330  cm  ~  '  for  the  HG2S 
was  only  13%  larger  than  the  corresponding  ratio  from  the 
T800.  This  value  is  just  outside  the  experimental  uncertainty 
of  ±  5%  in  the  base  line  determination  using  the  96% 
H,S04  comparison,  and  it  indicates  that  the  stray-light  level 
of  the  T800  is  slightly  better  than  that  of  the  HG2S,  as  noted 
previously.3 

The  total  integrated  BE  corrected  Raman  intensities,  /  / 


t  ’C 

FIG  2.  Integrated  Bove- Einstein  corrected  Raman  intensities  (O)  from 
water  vs  temperature  in  ’C. 


( 1  +  n),  i.e.,  the  areas  above  the  base  lines  in  the  region  0- 
325  cm  ~ 1  for  all  12  quantitative  spectra,  are  listed  in  Table  I. 
These  intensity  data  are  plotted  vs  temperature  in  *C  in  Fig. 
2,  where  the  linear  least  squares  fit  of  the  data  is  shown  by  a 
dashed  line.  The  linear  least  squares  equation,  not  shown  in 
the  figure,  is  /  =  22.636  —  0.070  847  t.  The  linear  least 
squares  slope  of  —  0.0708  47  represents  the  sole  significant 
feature  of  the  data.  This  slope  is  required  to  obtain  param¬ 
eters  relevant  to  the  YW  thermodynamic  method,33  dis¬ 
cussed  subsequently. 

B.  Thermodynamic  analysis 

1.  Determination  of  AH"  for  the  HB—NHB 
transformation 

Consider  the  unbonded  to  bonded  transformation 
U— B,  i.e.,  NHB— HB,  between  0-0  nearest-neighbors  in 
liquid  water,  where  U  refers  to  NHB  0-0  nearest  neighbors, 
e.g.,  severely  bent  and/or  stretched  O-H  O  interactions;  and 
B  refers  to  partially  covalent,  linear  and/or  weakly  bent, 
hydrogen  bonds,  O-H  —  O,  that  contribute  to  the  low-fre¬ 
quency,  0-350  cm-1,  Raman  intensity  from  water.  The  cor¬ 
responding  equilibrium  constant  is  [B]/[U],  However,  be¬ 
cause  of  conservation  of  all  0-0  nearest  neighbor 
concentrations,  [U]  +  [B]  =  a.  where  a  is  constant.34  The 
total  integrated  Raman  intensity  after  Bose-Einstein  correc¬ 
tion.  /,  is  considered  to  be  proportional  to  the  B  concentra¬ 
tion,  /  =  /?[B],  U  is  thought  not  to  produce  any  intermole- 
cu lar  Raman  contribution  whatsoever.  Substitution  of 
a  —  { B ]  for  [U  ].  and  1/(3  for  [B],  yields  an  equilibrium 
constant  of  the  form  I /(A  —  /],  where  A  =  ap.  The  tem¬ 
perature  dependence  of  the  equilibrium  constant  is  thus  giv¬ 
en 


where  A/f  °  =  -  RB,  and  AS"  =  RC. 

Young  and  Westerdahl  ( YW)  successfully  used  Eq.  ( 1 ) 
to  obtain  a  A H  °  value  from  Raman  data  for  a  two-state  equi¬ 
librium. "  Because ,0  cannot  be  measured  readily,  their  pro¬ 
cedure  w  as  to  choose  A  values  until  a  plot  of  In  [  /  /( A  -  />] 
yielded  linearity  in  l/T.  B  was  assumed  to  be  constant,  ie., 
A H  °  is  constant,  or  AC  P  equals  zero.  YW  then  determined  B 
and  C  from  least  squares  using  this  single  value  of  A.  The 
method  was  later  modified  by  varying  A  to  give  a  minimum 
in  the  least  squares  standard  error  of  fit.12  A  version  of  this 
modified  procedure  was  employed  initially  in  this  work  with 
the  data  of  Table  1,  and  it  yielded  A  values  between  26.4  and 
34.0.  However,  this  method  is  strongly  dependent  upon  the 
scatter  of  data  and  thus  a  new  procedure  was  adopted  here 
instead. 

Equation  ( 1 )  is  linear  in  terms  of  In  [7 /(A  -  /)  ]  vs  l/T, 
but  it  is  highly  nonlinear  in  terms  of  /  vs  T or  t.  From  the  data 
of  Table  1,  a  plot  of  Eq.  ( 1 )  using  least  squares  values  of  A,  B, 
and  C  (determined  by  methods  described  subsequently) 
yields  a  pronounced  sigmoid  shape  in  terms  of  I  vs  t.  Fig.  2. 
solid  curve.  The  maximum  value  of  the  absolute  slope  of  the 
sigmoid  curve  corresponds  to  the  inflection  point,  shown  by 
the  (perpendicular)  bar  at  about  41  *C.  This  slope  value 
gives  rise  to  the  new  method  of  determining  A ,  B,  and  C.  The 


essential  feature  of  the  method  involves  matching  the  slope 
at  the  inflection,  as  closely  as  possible,  to  the  linear  least 
squares  slope,  —  0.070  847,  from  /  =  22.636  —  0.070  847 1. 
See  the  juxtaposition  of  the  dashed  (linear  least  squares) 
straight  line,  and  the  solid  nearly  straight  segment  of  the 
sigmoid  curve  in  Fig.  2  between  about  —  25  and  100  'C. 

Slopes  of  sigmoid  curves  corresponding  to  Eq.  ( 1 )  plot¬ 
ted  in  terms  of  /  vs  t  are  given  by  the  following  equation 
which  results  from  differentiation  of  Eq.  ( 1 )  with  respect  to 
T: 

dl  dl  -  BHA  - 1) 

dt  dT  AT 2  1 

The  maximum  absolute  value  of  dl  /dt  corresponds  to  the 
slope  at  the  inflection  point.  Hence,  in  the  present  method, 
the  A  parameter  is  varied  until  the  slope  at  the  inflection 
approaches  the  linear  /  vs  /  least  squares  slope,  as  closely  as 
possible.  The  method  involves  finding  the  exact  inflection 
temperature,  and  then  the  corresponding  value  of  dl  /dt, 
from  Eq.  (2).  The  inflection  temperature  can  be  determined 
by  sol ving  Eq.  ( 1 )  for  and  Tm„  ( the  /  and  T  values  at  the 
inflection  point )  simultaneously  with  the  following  equation 
arrived  at  by  equating  d2I /dt 2  from  Eq.  (2)  to  zero, 

/  -  AT"'«  +  d_  (2) 

~  B  +  2  '  (3) 

Simultaneous  solutions  between  Eqs.  (1)  and  (3)  for 
and  are  obtained  by  iteration.  Then  7~lnfl  is  used  in  Eq. 

( 1 )  to  determine  the  slope  at  the  inflection,  which  is  the 
maximum  absolute  slope. 

From  least  square  treatment  of  the  Table  I  data,  and  the 
procedures  described,  a  plot  of  the  absolute  value  of  the  slope 
at  the  inflection  was  made  vs  A .  A  minimum  in  the  absolute 
slope  at  the  inflection  occurred  for  .4  =  28.15.  The  calculat¬ 
ed  Tin„  for  A  =  28. 15  is  4 1.25  °C,and/mn  is  20.82  °C,  see  Fig. 
2,  bar.  This  slope,  at  41.25  "C,  is  —  0.071  915,  which  only 
1.5^  larger  than  the  linear  I  vs  t  least  squares  slope  of 
—  0.070  847.  This  agreement  represents  the  best  possible 
agreement  between  the  two  different  least  squares  ap¬ 


proaches  to  the  Table  I  data,  because  no  inflection  slope  oc¬ 
curs  below  -  0.071  915  for  other  A  values  chosen.  More¬ 
over,  the  solid  curve,  and  the  dashed  line  of  Fig.  2  are 
visually  indistinguishable  between  about  —  25  and  100  *C, 
and  the  slopes  of  the  two  are  visually  identical  in  this  region. 

With  regard  to  linear  (/  =  22.636  —  0.070  847  /)  as  op¬ 
posed  to  nonlinear  [In (I /A  —  /)  =  B/T  +  C]  least  squares 
fits  fo  the  present  data,  it  should  be  emphasized  that  there  is 
little  difference  between  the  two  fits  for  the  temperature 
range  of  the  measurements,  Fig.  2.  But,  extrapolation  leads 
to  physically  unrealistic  results  for  the  linear  case,  as  op¬ 
posed  to  the  nonlinear  case.  For  example,  recent  Raman 
studies  of  shock-compressed  water  indicate  that  partially  co¬ 
valent,  charge  transfer  hydrogen  bonding  persists  to  1437  °C 
at  257  kbar,'5  whereas  the  linear  I  vs  t  equation  implies  that 
such  hydrogen  bonding  is  totally  absent  above  ~  320  °C.  The 
nonlinear  equation,  on  the  other  hand,  leads  to  the  (qualita¬ 
tively)  reasonable  value  of  0.09  for  the  mole  fraction  of  HB 
0-0  pairs  at  1437  °C.  Similarly,  the  linear  fit  requires  the 
amount  of  hydrogen  bonding  to  increase  linearly  below 
—  50  °C,  whereas  the  nonlinear  fit  indicates  that  hydrogen 
bonding  is  essentially  complete  near  50  °C.  In  this  regard,  it 
is  well  known  that  pure  water  cannot  be  supercooled  below 
about  —  40  *C.J6  This  fact  agrees  with  the  nonlinear  fit,  but 
not  with  the  linear  fit,  because  homogeneous  nucleation  is 
almost  certainly  related  to  the  absence  of  NHB  defects;  see 
the  sharp  bend  in  the  solid  curve.  Fig.  2,  below  —  40  °C. 
Thus,  the  importance  of  the  linear  fit  is  solely  related  to  its 
use  in  obtaining  parameters  for  the  nonlinear  YW  equation. 

The  A//°  for  the  reaction  U— •  B  obtained  from 
A=  28.15  is  —  2.6 ,  keal/mol  OHO,  or  —  2.6  ±  0. 1  kcal/ 
mol  OHO  when  all  of  the  experimental  uncertainties  are 
considered.  The  corresponding  A 5°  is  —  6.2,  cal/deg  mol 
OHO. 

Values  of  A//°  that  are  considered  to  be  reliable  and 
which  result  from  spectroscopic  measurements  are  listed  in 
Table  II.  The  quantitative  agreement  between  the  present 
A//°  value  of  2.6  +  0. 1  keal/mol  OHO  corresponding  to  the 
intermolecular  0-0  stretching  region,  with  A H  °  values  from 


TABLE  II  Values  of  enthalpy  (or  energy)  of  hydrogen-bond  rupture  for  waier  obtained  from  various  meth¬ 
ods 
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the  one-  and  two-phonon  OH-stretching  regions  clearly  in¬ 
dicates  that  the  same  process  is  involved.  The  present  A/7* 
value  is  tested  by  means  of  the  heat  capacity,  and  the  AS* 
value  is  examined  in  subsequent  comparisons. 

2.  Determination  of  extra  AH’ needed  for  complete 
separation  of  NHB  0-0  pairs 

The  extra  A H  *  required  for  complete  separation  of  NHB 
0-0  pairs  resulting  from  O-H  •••  O  breakage  is  readily  ob¬ 
tained  from  the  standard  enthalpy  change  for  vaporization 
at  0*C.  A H\  —  10.767  kcal/mol  H-O.45  This  value  corre¬ 
sponds  to  5.384  kcal/mol  0-0  pair,  because  there  are  two 
0-0  pairs  per  H-O  molecule  in  a  tetrahedral  network. 
( Here  no  distinction  need  be  made  between  HB  or  NHB  O- 
O  pairs. )  Thus, 

A//',  =/B(A//',  +  A//'-)  +  (1  -/e)A//'-,  (4) 

where  A H  \  is  the  enthalpy  of  vaporization  per  mole  0-0, 
A H\  is  2.6  kcal/mol  0-0,  and  A//,  is  the  extra  enthalpy 
required  for  the  complete  separation  of  the  NHB  0-0  pairs 
per  mole. 

The  value  of fB  at  0  °C  from  the  nonlinear  least  squares 
fit.  Fig.  2,  is  given  by  7(0  *0/28.15,  thus  fB  (0  °C)  =  0.842. 
Accordingly,  A H\  =  3.2  kcal/mol  0-0.  and  this  value  in¬ 
dicates  that  highly  cohesive  interactions  are  involved 
betw  een  the  O-H  and  the  O  in  these  0-0  pairs.  Note  that  the 
total  enthalpy  necessary  for  the  complete  separation  of  the 
intact  hydrogen  bonds,  i.e.,  to  break  the  bonds  and  then  to 
move  the  O-H  away  from  the  O  as  in  the  gas  at  0  *C,  is  5.8 
kcal/mol  0-0,  which  is  close  to  the  value  usually  quoted  for 
the  hydrogen  bond  enthalpy  or  energy.4''  However,  in  this 
work  it  is  emphasized  that  a  two-step  process  is  involved. 
The  first  step  is  defined  here  as  hydrogen  bond  rupture.  The 
second  process,  as  defined  here,  is  not  considered  to  involve  a 
hydrogen  bond  per  se,  as  seen  for  example  from  the  lack  of 
intermolecular  Raman  intensity.  A  more  precise  picture  of 
the  HB  and  NHB  classes  of  interaction  is  presented  subse¬ 
quently. 


3.  Calculation  of  Cp 

C  P  for  liquid  water  may  be  calculated  using  a  relation 
derived  by  Frank.47  According  to  this  derivation  the  stan¬ 
dard  heat  content  of  liquid  water  considered  as  an  ideal  solu¬ 
tion  is  given  by 

(5) 

or  by 

H)=Hu+fBUr.  (6) 

where  A//’  =  H B  —  H u  Differentiation  of  Eq.  (6)  with 
respect  to  T  yields  the  relation 

.  c,=c;>+/8ac;  +  a//*(^^.  (7) 

From  K  =  fB/{  1  —  fB  )  and  (3  In  K/3T)p  =  AT/  VRT\ 
one  obtains  (3ft/3T)p  =/,(  1  - f„)AH'/RT2.  Substitu¬ 
tion  for  ( 3fg/3T)p  in  Eq.  (7)  leads  to  the  result 


p  ~  C.  -*-/»  AC,  + 


/.(l  -fs) 


(“7- 


The  YW  method  described  previously  requires  A/7  *  to  be 
constant  over  the  temperature  range  involved,  and  thus 
A C\  =  0.  Hence,  the  second  term  in  Eq.  (8)  is  made  zero, 
for  consistency.  This  term  cannot  exceed,  and  is  thought  to 
be  smaller  than,  0.465  cal/deg  mol  0-0,  because 
C’p(s)  —C'p( g)  is  0.552  cal/deg  mol  0-0  at  0°C,  and 
fB  =  0.842.  The  third  term  of  Eq.  (8)  is  6.09  cal/deg  mol  O- 
O,  or  1 2. 1 8  cal/deg  mol  H-O.  The  standard  heat  capacity  of 
the  vapor  at  0”C  is  7.905  cal/deg  mol  H-O.  Hence  C\  is 
20.08  cal/deg  mol  H:0  or  1.11  cal/deg  g  H-O. 

The  calculated  heat  capacity  of  1.11  cal/deg  g  H-O  at 
0  °C  is  in  very  satisfactory  agreement  with  the  experimental 
value  of  1 .01  cal/deg  g  H-O,  when  it  is  considered  that  ideal 
solution  thermodynamics  was  employed,  and  more  impor¬ 
tantly  when  it  is  realized  that  errors  in  A/7*  are  squared  in 
Eq.  (8).  It  should  also  be  mentioned  that  the  specific  volume 
of  water  rises  3.6%  over  the  temperature  range  of  the  present 
experiments.  This  means  that  relatively  fewer  H:0  mole¬ 
cules  are  excited  in  the  path  of  the  laser  beam  after  expan¬ 
sion.  However,  a  density  correction  was  not  made  here  be¬ 
cause  of  the  scatter  of  the  data,  and  because  the  magnitude  of 
this  correction  is  about  the  same  as  the  long  term  drift  of  the 
laser  power.  Nevertheless,  this  correction  would  be  in  the 
direction  of  lowering  the  A/7°  value,  which  would  improve 
the  calculated  C"  value. 

Pi 

4.  Entropies 

Consider  that  S\  for  liquid  water  is  given  by  the  ideal 
solution  formulation 

S  I  ~  ZbS  B  +fuSu  —  f B  1 nfB  -r  fu  Inf,),  (9) 
where  the  last  term  is  the  entropy  of  mixing.  Substitution  of 
A S’  =  S B  —  S u  in  Eq.  (9)  yields 

S]  =S'B  +  (fB-  DA S‘-Jt(fB  lnfB  +f  In/.). 

(10) 

At  0°C,/g  =  0.842  and  fu  =0.158.  AS*  =  -  6.21  cal/deg 
mol  0-0.  The  accepted  value  of  S ,  at  0  °C  is  15.166  cal/ 
deg  mol  H-O,41*  or  7.583  cal/deg  mol  0-0  pair.  Hence,  from 
Eq.  ( 10),  SB  =  1 1  47  cal/deg  mol  H-O.  But  the  standard 
entropy  of  ice  at  0  *C  is  9.909  cal/deg  mol  H-O.  Therefore, 
S B  is  larger  by  1.56  cai/deg  mol  H-O.  This  —  16%  addi¬ 
tional  entropy  means  that  the  bonded  component  of  water  is 
more  disordered  than  ice.  and  thus  cannot  be  identical  to  ice. 

One  possible  source  of  extra  disorder  of  the  HB  compo¬ 
nent  may  be  that  some  or  most  of  the  hydrogen  bonds  in  the 
bonded  form  of  water  are  not  quite  linear,  i.e.,  a  range  of  HB 
O-H  —  O  angles  somewhat  less  than  180°  exists. 

The  standard  entropy  for  the  NHB  form  of  liquid  w  ater 
is  also  readily  calculated  at  0*C.  —A S°  =  SU~SB 
=  12.42*  cal/deg  mol  H-O.  S B  =  1 1.47  cal/deg  mol  H;0. 
Hence  S  u  =  23.90  cal/deg  mol  H-O.  How  ever,  the  spectro¬ 
scopic  value  of  S*  for  the  gas  at  25  *C  is  45.106  cal/deg  mol 
H-O,4"  which  upon  adding  Cp  In ( 273. 15/298. 15)  yields  an 
S°  value  of  —44.4  cal/deg  mol  H-O  at  0  ’C.  Because  this 
value  is  larger  than  S'u  by  -20.5  cal/deg  mol  H-O,  it  is 
obvious  that  the  NHB  component  is  much  more  ordered  and 
much  more  liquid-like  than  the  gas.  This  order  is  related  to 
the  strong  cohesive  forces  which  refer  to  the  AT/',  value  of 


(8) 


3.2  kcal/rrioi  0-0  pair,  and  probably  arises  from  strong 
electrostatic  ( multipole ) ,  induction,  and  dispersion  forces.41' 

Finally,  it  should  be  noted,  that  although  S  *  ( g,  273. 1 5 
K)  >5 it  is  nevertheless  reasonable  to  use  the  gas  value  for 
C'  mEq.  ( 8 ).  The  heat  capacity  of  ice  at  0°C  is  nearly  equal 
to  the  value  for  the  vapor  at  that  temperature,  i.e.,  only  high¬ 
er  bv  1.104  cal/deg  mol  H:0.  The  dominant  factor  in  the 
high  heat  capacity  of  the  liquid,  compared  either  to  ice  or  to 
steam,  is  clearly  the  presence  of  the  U— B  or  NHB— »HB 
equilibrium  and  the  large  A//’,  associated  with  it. 

C.  Contour  analysis  and  mode  softening 

Computer  decompositions  of  the  low -frequency  con¬ 
tours  were  accomplished  by  use  of  an  analog  computer.  Ra¬ 
man  amplitudes  ( the  differc ce  between  the  total  amplitude 
and  the  base  line)  were  transferred  manually  to  horizontal 
base  lines  using  frequency  increments  of  12  5  cm"  Analog 
decomposition  was  then  earned  out  using  two  (nominal  60 
and  170  cm  ')  Gaussian  components.  It  was  found  in  all 
cases,  however,  that  the  component  intensity  below  60  cm  ~  1 
had  to  be  truncated  such  .hat  it  not  exceed  the  observed 
intensity  and  so  that  its  tail  not  be  allowed  to  extend  (slight¬ 
ly  )  into  the  ann-Stok  .s  region.  The  60  cm  ~ 1  component  that 
resulted  was  thus  o\iy  Gaussian  in  shape  above  its  peak,  i.e., 
it  could  be  considered  to  be  skewed  to  high  frequencies.  Also 
neutron  inelastic  scattering  data  y  ield  a  peak  in  the  spectral 
density  function  of  water  at  —  100  cm  -  1  w  ith  a  shoulder  at 
60  cm  1  This  neutron  observation  might  explain  thenon- 
symmetnc  Raman  shape  seen  here  for  the  60  cm  ~  1  compo¬ 
nent  Furthermore,  a  nearly  identical  nonsymmetric  Raman 
shape  has  been  seen  for  an  isolated  low-frequency  compo¬ 
nent  of  molten  B.O,  at  —28  cm'  With  the  two  compo¬ 
nents  described,  a  consistent  sei  of  component  peak  frequen¬ 
cies  and  full  widths  at  half-height  ( FW'HH )  was  obtained 

The  component  peak  frequencies  are  plotted  vs  tem¬ 
perature  in  °C  in  Fie  3  The  FWHH  values  (not  shown) 
were  nearly  constant.  —65  cm  '  for  the  60  cm  1  compo¬ 
nent.  and  —  140  cm  1  for  the  170  cm  1  component 

Decreases  in  the  component  frequencies  are  evident 
from  Fig.  3.  Least  squares  equations  referring  to  these  de¬ 
creases  are  shown  in  the  figure  The  frequency  change  of  the 
170  cm  1  component  agrees  satisfactorily  with  recent  re¬ 
sults.-' 

The  frequency  decrease  observed  for  the  1 70cm  '  com¬ 
ponent  with  temperature  rise  is  consistent  with  a  decrease  in 
the  stretching  force  constant  of  the  partially  covalent  FI  -  O 
bond  of  the  FIBO-H  —  O  unti.  Such  a  decrease  in  the  FI  —  O 
stretching  force  constant  is  consistent  with  an  increase  in  the 
H  •••  O  bond  distance,  as  seen  from  the  inverse  relationship 
between  bond  distance  and  force  constant. ' 1  Flow  ever,  the  x- 
ray  RDF  shows  little  increase  the  nearest-neighbor  0-0 
distance  over  the  small  temperature  range  involved  here.': 
Hence,  an  increase  in  the  H  •••  O  distance  would  seem  to 
demand  a  decrease  in  the  O-H  •••  O  angle,  if  the  0-0  dis¬ 
tance  remains  essentially  constant. 

The  frequency  decrease  observed  for  the  60  cm  ' 1  com¬ 
ponent  also  indicates  that  a  decrease  in  some  force  constant 
has  occurred,  probably  the  O-H  O  bending  force  con¬ 
stant  An  increase  in  the  O-H  —  O  angle  in  the  direction  of 
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FIG  3  Raman  comri'nem  frequencies  ( two  component  analysis)  vs  tem¬ 
perature  in  "C.  from  Bose- Einstein  .orrected  spectra  for  water 


180°  would  thus  normally  be  expected,  because  this  effect 
would  lower  the  0-0  repulsion  by  moving  the  oxygen  atoms 
farther  apait.  cf.  Refs.  53.  54  for  vitreous  SiO:.  But  an  in¬ 
crease  in  the  O-H  O  angle  is  extremely  unlikely  ,  because 
more  nearly  linear  hydrogen  bonds  are  generally  thought  to 
he  stronger,  not  weaker  Moreover,  this  explanation  for  the 
60  cm  '  frequency  decrease  is  just  the  opposite  of  that  for 
the  170  cm  1  frequency  decrease.  Obviously,  another  more 
fundamental  effect  is  involved. 

It  is  v  lriually  certain  that  the  occurrence  of  intermolecu- 
lar  Raman  intensity  for  water  is  the  result  of  the  partial  cova¬ 
lency  of.  or  of  the  charge  transfer  into,  the  H  —  O  part  of  the 
O-H  •••  O  units.4 '  "  Furthermore,  this  partial  covalency  has 
been  show  n  to  decrease  as  the  O-H  O  angle  decreases  be¬ 
low  180°.*  although  a  concomitant  increase  in  the  H  •••  O 
distance  would  also  have  similar  effects.  At  any  rate,  a  de¬ 
crease  in  the  partial  covalencv  of  the  H  —  O  bond,  would 
decrease  both  the  O-H  —  O  stretching  and  bending  force 
constants  seen  here,  because  both  the  60  and  the  170  cm  - 1 
components  are  Raman  active  by  virtue  of  this  partial  cova¬ 
lency.  Furthermore,  when  the  O-H  •••  O  angle  becomes  very 
small,  say  1 50°  or  less,  and  when  the  H  —  O  distance  is  large, 
for  example,  -2.4  A,  no  intermolecular  Raman  intensity 
whatsoever  would  be  expected.  Sec.  Ill  D.  In  this  case  other 
types  of  cohesive  interactions  would  be  involved.  Secs. 
Ill  B  2,  III  B  4.  and  III  D,  where  an  NHB  OH-stretchmg 
frequency  near  3650  cm  ~  1  results. 


INTENSITY 


Three  component  analyses  of  the  low-frequency  con¬ 
tour  were  also  earned  out  using  the  truncated  or  skewed  60 
cm~'  component,  plus  two  other  completely  symmetric 
Gaussian  components.55  The  three  component  fits  were 
somewhat  better  than  those  described,  but  they  do  not  alter 
the  preceding  conclusions. 

A  consequence  of  the  O-H  •••  O  bending  hypothesis  is 
that  the  present  thermodynamic  treatments  must  involve 
average  A/f *  and  AS*  values,  where  these  averages  refer  to 
the  HB  and  NHB  population  groups.  Nevertheless,  it  is  evi¬ 
dent  that  consistent  definitions  of  these  groups  are  involved, 
because  of  the  quantitative  agreements  between  A H°  values 
from  the  tntermolecular  and  the  one-  and  two-phonon  OH- 
stretching  regions  from  both  Raman  and  infrared  measure¬ 
ments.  as  well  as  other  nonvibrational  types  of  measure¬ 
ments,  Table  II. 

The  average  nature  of  the  and  AS  ’values  in  terms  of 

OH-stretching  vibrations  refers  to  the  large  spectral  widths 
observed.  The  two  population  classes  are  evident  from  the 
loss  of  Raman  intensity  in  the  HB  OH-stretching  region  and 
the  concommitant  intensity  gain  in  the  NHB  OH-stretching 
region  with  temperature  rise.  The  regions  are  separated  by 
an  isosbestic  point  which  occurs  when  the  compound  fre¬ 
quencies.  and  thus,  the  corresponding  average  thermody¬ 
namic  quantities,  do  not  change  rapidly  with  temperature. 
These  spectroscopic  features  of  water  are  all  involved  in  the 
quantitative  X(Z.Z  +  A-)  Y  OH-stretching  and  in  the  corre¬ 
sponding  difference  spectra  described  next. 

D.  Isosbestic  point  and  Raman  difference  spectra  from 
the  fundamental  OH-stretching  region 

Raman  spectra  obtained  between  3.5  and  72.0  C  with 
the  HG2S  double  monochromator  under  precisely  constant 
conditions  of  excitation  and  detection  are  shown  in  Fig.  4. 
The  laser  beam  position  in  the  filled  Raman  cell  was  held 
exactly  fixed  at  all  temperatures  by  use  of  New  ton's  rings,  as 
described  in  Sec.  II.  The  spectra  shown  in  Fig.  4  were  ob¬ 
tained  with  dc  detection  and  refer  to  the  X(  ZJl  —  Z)  )'  ge¬ 
ometry.  which  was  used  throughout  this  work.  [Quantita¬ 
tive  Raman  spectra  corresponding  to  the  ,V(  YJC  -r  Z)Y 
X(ZZi  Y.  and  XiZX)  Y  geometries  for  water  are  presented 
in  the  preceding  article.5'’] 


The  X(ZJC  +  Z)Y  spectra  of  Fig.  4  display  a  precise 
isosbestic  point  at  3425  cm-1.  No  evidence  whatsoever  is 
present  in  Fig.  4,  or  in  the  large  original  spectral  tracings  that 
would  indicate  that  the  crossings  occur  within  a  sizable  area; 
or,  that  the  crossings  occur  at  the  same  frequency,  but  not  at 
the  same  intensity.  [  Precise  isosbestic  points  were  also  ob¬ 
served  for  the  X(  YX  -f-  Z)  Y,  X(ZZ)  Y,  and  X{ZX)  Y  spec¬ 
tra,  and  a  precise  isosbestic  point  was  observed  for  the  a 2 
spectrum  as  well,  e.g.,  from  /(VV)  —  4/3/ ( VH).56  This 
precise  new  a 2  isosbestic  point  indicates  that  the  data  of  Ref. 
32  are  in  error,  j 

The  isosbestic  point  of  Fig.  4  may  be  understood  by  the 
following  analysis  which  employs  HB  and  HHB  interac¬ 
tions. 

Consider  that  the  Raman  intensity  l{co,t),  where  co  re¬ 
fers  to  frequency  and  t  to  temperature,  is  given  by 

/  (fij,/)  =  /HB(r)-t-/NHB(r),  (11) 

_  ,  ,,-._-x.  ■  t 

or 

/(w)  =  2^Q(7)  +  I^;c;(  o,  (i2) 


where  the  index  cj  can  run  over  the  entire  OH-stretching 
contour  to  allow  for  component  overlap,  J u  is  a  constant  for 
a  particular  aj  and  refers  to  HB  interactions,  and  J'tj  is  the 
corresponding  constant,  but  refers  to  NHB  interactions,  and 
C„  and  Ci  are  temperature  dependent  HB  and  NHB  con¬ 
centrations,  respectively.  Also 

veto  +  vc:(f)  -  a,  (13) 

<±>  ui 

where  a  is  a  constant.'4  i.e..  conservation  of  mass  applies.  At 
the  isosbestic  point,  co  —  /,  I{co,t)  =  const,  and  dI(co,i)/ 
di  =  0.  Hence, 


,  dC At) 

J, - j  ; 

dt 


1£2±  =  o. 

dt 


and 


dC,(t )  dC’tn 

dt  dt 


(14) 

(15) 
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KIG  4  Quantitative  ( absolute ) 
Raman  spectra  in  the  OH 
stretching  region  obtained  under 
rigorously  constant  conditions  of 
excitation  iniensity  and  detec¬ 
tion  geometry  Note  the  isosbes¬ 
tic  point  at  3425  cm  " '.  where  all 
spectra  have  the  same  intensity 
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FIG  5  Raman  difference  spectra  from  water  in  the  fundamental  (one- 
phonon)  OH-stretching  region.  Note  the  isosbestic  frequency  at  3425 

cm-1 


Substitution  ofEq.  (14)  into  Eq.  (15)  yields 
dC,(t) 

—  =  o.  (16) 

dt 

But.  dC  ( t)/dt  cannot  reasonably  be  zero/  Therefore. 
J,  =  J at  the  isosbestic  point.  As  both  concentrations  C,(r ) 
and  C'(t)  change  with  temperature,  the  equality  of  J ,  and 
J '  requires  the  Raman  intensity  at  the  isosbestic  point  to 
remain  constant,  because  the  intensity  per  mole  liter  is  the 
same  for  either  species  or  interaction.  Furthermore,  it 
should  be  made  clear  that  the  equality  of  J,  and  J '  definitely 
does  not  mean  that  the  two  species  or  two  classes  of  interac¬ 
tions  refer  to  the  same  structure.  On  the  contrary,  the  inter¬ 
action  energy  between  H-O  molecules  is  a  function  of  the  O- 
O  separation,  of  two  polar  angles,  and  three  Euler  angles.'8 
This  means  that  the  NHB  and  HB  classes  of  interaction  may 
be  regarded  as  separate  volume  elements  in  this  6  space. 
Hence,  in  the  specific  case  of  the  isosbestic  point,  mapping  of 
the  6  space  into  the  intensity  vs  frequency  space  leads  to  a 
situation  in  which  two  different  interactions.  HB  and  NHB, 
have  the  same  molar  intensity,  J  =  J at  oj  = 

Raman  difference  spectra  were  obtained  with  the  T800 
triple  monochromator  for  temperatures  from  3.5  to  89.3  ”C. 
(The  noise  levels  are  larger  than  those  of  Fig.  4  because  the 
conditions  were  different,  and  amplification  of  the  differ¬ 
ences  was  also  involved.)  The  intensity  differences  are  all 
zero  at  3425  cm  ' ',  as  shown,  for  example,  by  spectra  corre¬ 
sponding  to  three  temperatures  (of  a  total  of  12)  in  Fig.  5. 
The  difference  spectra  of  Fig.  5  were  obtained  using  digitized 
data  and  the  computer  described  in  Sec.  II  From  Fig  5  it  is 
evident  that  the  integrated  intensity  loss  below  3425  cm  '. 
upon  temperature  rise,  is  accompanied  by  a  somewhat 
smaller  intensity  gain  in  the  integrated  intensity  above  3425 
cm  '  .  This  observation  suggests  that  (he  HB  OH-stretching 
oscillators  are  replaced  by  NHB  OH-stretching  oscillators  at 
a  higher  frequency.  Moreover,  the  molar  intensity  from  the 
NHB  OH-stretching  oscillators  must  be  relatively  smaller 
than  those  of  the  HB  oscillators  But  before  it  can  be  con¬ 
cluded  that  there  is  an  exact  conservation  between  oscilla- 
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FIG  6  Integrated  Raman  intensity  gain  above  the  isosbestic  point  vs  tnte- 
grated  intensity  loss  below,  see  Fig.  5.  X(Z.Z  +  X)  >' 


tors  of  the  two  classes,  it  is  necessary  to  examine  the  intensity 
differences  quantitatively. 

The  NHB  integrated  intensity  gain  above  3425  cm  ~  1  is 
plotted  vs  the  HB  integrated  intensity  loss  below  3425  cm  “ 1 
(absolute  values)  in  Fig.  6.  Proportionality  is  definitely  evi¬ 
dent  from  Fig.  6  for  the  temperature  range  of  3.5  to  89.3  °C 
The  linear  least  squares  fit  of  the  data  ( including  the  point  at 
the  origin )  is  shown  by  a  straight  line. 

The  proportionality  evident  from  Fig.  6  may  be  seen  to 
arise  from  conserv  ation  of  HB  and  NHB  OH-stretching  os¬ 
cillators.  Consider  then  the  integrated  Raman  intensity  be¬ 
low  3425  cm  “ 1  in  the  OH-stretchmg  region,  /HB  ( t ) ,  is  given 
bv 


,  U)  =  =yHB  V  C„(r), 


(17) 


where  C„,  (O  is  the  temperature  dependent  concentration  of 
3H-stretching  oscillators  at  a  single  frequency  w.  and7HB  is 
an  average  scattering  per  total  mole  of  OH  oscillators  per 
iter  as  defined  by  the  second  equality  of  Eq  (17).  The  inte¬ 
grated  Raman  intensity  above  3425  cm  -  1  in  the  NHB  OH- 
vtretching  region.  /NHB  (r),  is  similarly  given  by 


,(t)  =  ^y.:,c:(r)  =yN.HB  vc;,(f). 


(18) 


Substitution  of  Eqs  (17)  and  (18)  into  Eq.  (13),  which 
involves  conservation  of  the  total  HBand  NHB  OH-stretch- 
ing  oscillator  concentrations,  and  rearrangement  yields 


,(/) 


-  (  J"n  V 

~\  7  r  NMB 

Wshh  ' 


tO-  aJH 


(  19) 


A  difference  equation  results  from  Eq.  (19)  for  tempera¬ 
tures  n  and  r„  where  r,  =  const  =  (y),  thus 

-  [W':>  ~  WD] 

=^— - 1  [fsHB  (  )  _  ^NHB  (  Y)  ]  (20) 

•^SHB  I 

From  Eq.  (20)  it  is  obvious  that  the  absolute  value  of  the 
integrated  difference  -  [/„B(r:)  -  increases  in 

direct  proportion  to  the  increase  in  the  other  integrated  dif¬ 
ference.  [/nhbO:)  _  ^nhb  (y)  ] *  as  t  increases  above  y. 

In  regard  to  Eq.  (13)  and  to  the  [  U  ]  +  [  B  ]  =  a  conser¬ 
vation  condition  in  the  YW  derivation,  it  should  be  made 
clear  that  because  the  A H°  value  from  the  OH-stretching 
region  agrees  with  the  low-frequency  intermolecular  A H° 
value,  it  is  apparent  that  the  same  two  population  classes  HB 
and  N'HB.  must  be  involved  in  both  cases.  Hence.  Eq.  (13)  is 
simply  a  more  detailed  expression  of  [U]  +  [B]  =  a.  i.e.. 

~  [ U ]  and  C„  ( t )  =  [ B ) .  Moreover,  this  iden¬ 
tity  of  definitions  must  be  demanded,  if  exactly  the  same 
process  is  involved  in  the  intermolecular  and  intramolecular 
Raman  measurements.  A  further  consequence  of  the  formu¬ 
lation  of  Eq.  (13)  is  that  a  range  of  structures  must  exist  for 
HB  and  NHB  population  classes.  However,  this  range  of 
structures  in  no  way  precludes  the  present  thermodynamic 
analysis,  because  the  A H°  and  A S‘  values  refer  to  averages 
over  the  configurations  involved. 

It  should  also  be  noted  that  the  conclusions  resulting 
from  Figs.  4,  5,  and  6,  and  in  particular  Eq.  (20),  provide 
definite  e\  idence  for  the  presence  of  the  HB— NHB  transfor¬ 
mation.  but  they  fail  to  provide  any  evidence  whatsoever  for 
Fermi  resonance,5"'  or  for  untested  and  contrived  explana¬ 
tions.  such  as  electrical  anharmomcity.5’0-'’ 1  Futhermore,  the 
present  data  clearly  eliminate  assignments  of  the  3630  cm  ~ 1 
shoulder  in  the  OH-stretching  region  to  a  combination  vi¬ 
bration.''4  On  the  other  hand,  the  data  of  Figs.  4,  5,  and  6 
agree  completely  with  molecular  dynamics,  cf.  Fig.  1 1  of 
Ref.  63. 

E.  Estimation  of  NHB  O-H  O  angles  and  distances  using 
vibrations  from  vicinal  surface  silanol  groups  on 
vitreous  silica 

The  Raman  spectrum  from  dilute  HDO  in  D  O  ( I 
mol  H  O  in  D  O.  25  °C)  displays  a  broad  shoulder  compo¬ 
nent  centered  ner  3628  cm*  1  due  to  NHB  OH  stretching,  in 
addition  to  the  mam  peak  at  —3435  cm  1  from  HB  OH 
stretching.5’5"  This  HDO  spectrum  is  useful  in  the  following 
considerations. 

Silanol  groups,  Si-O-H  ( three  Si  bonds  not  expressed  ). 
like  HDO  and  HTO  molecules. 51  involve  decoupling  of  the 
stretching  vibrations.  The  Si-O-H  decoupling,  however,  is 
more  complete  On  a  vitreous  silica  surface,  the  free  ( dis¬ 
play  s  P  and  R  branches5’*)  OH-stretching  frequency  of  Si- 
O-H  is  3750  cm  ~  ” and  the  Si  vs  OH  stretching  frequen¬ 

cy  is  970-980  cm  171  However,  when  the  surface  Si-O-H 
coverage  becomes  large,  a  broad  peak  develops  near  -  3650 
cm’1,  in  addition  to  the  sharp  3750  cm “  1  peak.5'"  "  This 
3650  cm  1  peak  has  been  shown  conclusively  to  arise  from 
NHB  or  O-H  O  interactions  between  neighboring  or  vicinal 


silanol  groups  whose  Si  atoms  are  linked  by  a  surface  Si-O- 
Si  bridge.'’*-70 

From  known  angles  and  bond  distances  it  is  possible  to 
estimate  the  O-H  O  angle  and  H  O  distance  corresponding 
to  the  3650  cm  ~ 1  feature.  For  bulk  fused  silica77  or  for  thin 
surface  silica  films75  the  Si-O,  O-O,  and  Si-Si  nearest-neigh¬ 
bor  distances  are  1.6,  2.6,  and  3. 1  A,  respectively,  and  the  Si- 
O-Si  bridging  angle  is  144°. 74  75  The  tilt  or  dihedral  angle  <5 
between  bridged  Si04  tetrahedra  is  about  16°, 75  and  the  Si- 
O-H  angle  /?  is  near  1 13°  (experiment)5’*  or  about  109°  the¬ 
ory,75’  i.e.,  both  values  are  not  far  from  the  tetrahedral  angle. 
109°28'.  The  O-H  distance  is  0.96  A.7"’  From  these  data,  the 
calculated  O-H  O  angles  and  H  O  distances  are,  respective¬ 
ly.  148°  and  -2.4  A  (0  =  113°,  <5=  16°),  and  153°  and 
—  2.4  A  (  (5  =  109,  6  =  16°).  The  calculated  secon4-neigh- 
bor  O'-O'  distance  (O'-Si-O-Si-O')  is  3.2  A,  but  this  is 
actually  the  0-0  distance  of  the  O-H  O  interaction  involved 
here. 

The  nearest-neighbor  0-0  distance  for  ice  Ih  is  2.76 
A,’*  and  theH  —  Odistanceis  thus  about  1.8  A  when  theO- 
H  bond  distance  of  0.96  A  is  employed.  For  comparison,  the 
H  O  distance  corresponding  to  the  3650  cm  ~ 1  vicinal  silanol 
frequency  is  almost  as  large  as  the  entire  nearest-neighbor 
0-0  distance  in  ice  Ih,  that  is.  the  vicinal  H  O  distance  is 
33%  larger  than  the  ice  Ih  H  —  O  distance.  Furthermore, 
the  vicinal  O-H  O  angle  has  changed  by  30°  compared  to  the 
ice  Ih  angle;  a  large  decrease  for  H-bonded  systems.  Accord¬ 
ingly.  the  vicinal  O-H  O  interaction  can  hardly  be  consid¬ 
ered  as  a  hydrogen  bond.79  and  the  1 50°  O-H  O  angle  and  the 
2.4  A  H  O  distance  must  certainly  refer  to  a  weakened  inter¬ 
action,  clearly  falling  w  ithin  the  present  definition  of  nonhy¬ 
drogen-bonded  interactions.  Such  large  angle  and  distance 
deviations  from  the  ice  lh  values  would  be  expectd  to  de¬ 
crease  the  partially  covalent  or  charge  transfer  character  of 
the  H  —  O  linkage. J_h  and  would  also  be  expected  to  lower 
the  OH-stretching  intensity,  and  preclude  significant  inter¬ 
molecular  Raman  intensity  in  the  170  cm-  1  region. 

With  regard  to  relative  HB  and  NHB  OH-oscillator 
concentrations  it  should  be  noted  that  sizable  Raman  inten¬ 
sity  occurs  at  3650  cm  ~  1  for  dilute  HDO  in  D-O.5"5'  The 
contour  height  at  3650  cm  ’  is  -  1/4  the  peak  height  at 
3435  cm  " l.5'5’  Thus,  from  the  silanol  calculation  one  would 
expect  that  a  significant  fraction  of  the  interactions  in  this 
HDO-D.o  mixture  would  refer  to  O-H  O  angles  of  1 50°  and 
H  O  distances  of  2  4  A  Moreover,  a  similar  situation  would 
be  expected  for  liquid  water,  because  the  water  spectra  are 
not  changed  very  greatly  by  coupling  effects 

It  should  be  emphasized,  however,  that  the  entire  spec¬ 
tral  region  from  the  isosbestic  frequency  of  3425  cm  "  1  to  the 
limit  of  about  3750-3800  cm-1  refers  to  a  range  of  NHB 
O-H  O  angles  and  NHB  0-0  distances,  of  which,  the  calcu¬ 
lated  silanol  values  represent  a  single  contribution,  i.e.,  one 
complexion  of  the  ensemble 

Finally.  Narten  and  Levy57  have  made  accurate  x-ray 
measurements  of  G(r)  for  liquid  water  between  4  and 
200  °C  Their  <7(r)  data  for  4  °C  show  asymmetry  at  —  3.2- 
3.3  A  above  the  intense  nearest-neighbor  0-0  peak  at  2  9  A 
This  asymmetry  contribution  at  3  2-3  3  A  is  consistent  with 
the  present  conclusions  than  an  0-0  distance  of  -  3  2  A 


gives  rise  to  a  frequency  just  above  the  centroid  of  the  NHB 
OH-stretching  component,  which  is  3628  cm  ~ 1  for  HDO.  j 
The  temperature  dependence  of  the  x-ray  data  is  also  consis¬ 
tent  with  the  observation  that  the  relative  NHB  OH-stretch¬ 
ing  intensity  from  HDO  increases  with  temperature  rise,67 
i.e.,  that  the  concentrations  of  OH  oscillators  involved  in 
long  0-0  distnaces,  ~  3.2  A,  increases  with  increasing  tem- 
perture. 

IV.  SUMMARY 

BF  correction  of  low-frequency  Raman  spectra  from 
water  and  thermodynamic  analysis  of  the  resulting  integrat¬ 
ed  intensities  yields  a  A H  °  value  of  2.6  ±  0. 1  kcal/mol  OHO 
for  the  rupture  of  hydrogen  bonds.  This  AH  °  leads  to  reason¬ 
able  values  for  the  heat  capacity  and  for  various  entropies, 
and  is  in  agreement  with  AH°  values  obtained  from  vibra¬ 
tional  spectroscopy  for  the  one-  and  two-phonon  OH- 
stretching  regions,  as  well  as  from  other  methods  such  as 
depolarized  Rayleigh  scattering  or  molecular  dynamic 
Some  implications  of  the  new  result  are:  that  the  V  -- 
quency  mtermolecular  Raman  spectrum  arises  predomi¬ 
nantly  from  partially  covalent  linear  and/or  weakly  bent  O- 
H  —  O  hydrogen  bonds,  HB;  that  the  breakdown  of  such 
partially  covalent  interactions  is  responsible  for  the  high 
heat  capacity,  compared  to  ice  or  to  steam;  and  that  this 
breakdown  produces  another  class  of  cohesive  interactions, 
NHB.  whose  further  complete  separation  would  require  an 
additional  3.2  kcal/mol  0-0  pair. 

Quantitative  (absolute)  Raman  spectra  obtained  in  the 
OH-stretching  region  yield  a  precise  isosbestic  point  at  3425 
cm  '  '.  Such  a  point  is  indicative  of  an  equilibrium  between 
the  HB  and  NHB  classes  of  interaction.  Raman  difference 
spectra  in  the  OH-stretching  region  show  a  rise  in  the  total 
integrated  intensity  above  the  isosbestic  frequency,  and  a 
decrease  below.  A  plot  of  the  coresponding  integrals  indi¬ 
cates  direct  proportionality.  This  fact  establishes  the  pres¬ 
ence  of  the  HB— NHB  equilibrium,  and  indicates  that  the 
OH-oscillator  concentration  is  conserved  in  terms  of  the  HB 
and  NHB  classes. 

V.  ADDENDUM 

A.  Volume  changes  and  pressure  effects 

The  equation 
Inf/ /(28.15  -  /)) 

=  ln[/,/(l  — /*)] 

=  13.1 1 1  024  97"  -  3.126  419  25, 

see  Fig.  2,  allows  values  of  /  or  fB  to  be  calculated.  Values  of 
f8  are  useful  in  making  rough  estimates  estimates  of  the  vol¬ 
ume  change,  A  V  —  VB  —  F„ ,  corresponding  to  the  transfor¬ 
mation  U— B. 

The  volume  change  A  V  may  be  approximated  from  the 
ideal  solution  formulation 

K  =/, Vg  +/,  v„  =/,  f,  +  <  i  -/„ ) f„,  (21) 

which  in  terms  of  the  volume  change  A  F  is 

K  =  /» A F -r-  Vu.  (22) 


TABLE  III.  Values  of  the  partially  covalent  hydrogen-bonded.  O-H  —  O 
mole  fraction,/,,  at  various  temperatures  for  water. 


t.’C 

/. 

0 

0.842 

10 

0818 

20 

0  794 

25 

0  781 

30 

0.768 

40 

0.743 

50 

0.717 

60 

0  692 

70 

0  667 

80 

0  643 

90 

0619 

too 

0  596 

If  values  of fB  and  V ,  are  known  for  widely  spaced  tempera- 
simultaneous  solutions  of  Eqs.  (21 )  or  (22)  yield  AF. 

.  or  temperatures  of  4  and  90  °C,  which  span  the  present  low  - 
frequency  Raman  data,  simultaneous  solution  yields  an 
average  AF  value  of  —  1.5  cmVmol  OHO.  Alternatively, 
Eq.  (22)  may  be  differentiated  with  respect  to  T  or  to  P. 
Then  from  the  forced  assumption  that  all  volumes  except  F 
are  independent  of  T  or  of  P,  and  from  substitution  of  the 
appropriate  expression  for  (dfB/3T)Por  (3fB/3P)T,  seethe 
derivation  in  Sec.  Ill  B  3,  values  of  AFvs. /may  be  calculat¬ 
ed  from  known  values  of  (3V,/3T)P  or  (3V,/3P)T.  When 
( 3V,/3T)r  is  equated  to  A  FA HfBfu/RT',  the  values  of  AF 
that  result  are  as  follows:  ( 1 )  0  °C,  -r  0.26  cm  ‘/mol  OHO; 
( 2 )  20 °C,  -  0.75;  (3)  40 °C,  -1.4:  (4) 60°C,  -1.9;  and 
(5)  80  °C.  —  2.5.  Of  course,  because  these  A  Fvalues  change 
both  in  sign  and  in  magnitude  from  0  to  80  °C,  it  is  evident 
that  the  assumption  relative  to  temperature  independence  of 
AF,  VB,  or  Vu  is  faulty,  but  all  of  the  calculated  A  Fvalues 
are  at  least  small.  Moreover,  the  0  °C  result  agrees  with  a 
value  reported  by  Angeil.  Ref.  36.  The  important  point  here 
is  that  if  AF  is  small,  its  magnitude  is  in  qualitative  agree¬ 
ment  with  Ref.  80.  where  no  significant  change  in  the  NHB 
OD-stretching  component  intensity  relative  to  the  HB  OD- 
stretchmg  component  intensity  could  be  ascertained,  within 
the  accuracy  of  the  computer  analysis  employed,  to  roughly 
lOkbarat  32  °C. 

For  convenience  values  fo/fi  in  10°  increments  between 
0  and  100  °C  are  listed  in  Table  III. 

It  has  recently  been  possible  to  obtain  Raman  spectra  in 
the  fundamental  OH-stretching  region  from  strongly- 
shocked  water  at  pressures  of  257  kbar  and  temperatures  to 
1437  °C  ( 1710  K),  where  the  density  is  2  gem-3.’5  HB  and 
NHB  classes  or  interactions  are  both  definitely  indicated  up 
to  1437  °C  in  agreement  with  the  present  work.  However, 
most  of  the  hydrogen  bonds  are  broken  at  1437  eC,  and  the 
H.O  molecules  are  thought  to  be  held  together  by  other 
strong  forces,  e.g.,  electrostatic  (multipole),  induction,  and 
dispersion.  At  a  density  of  2  g  cm-3,  the  OH-stretching  vi¬ 
brations  of  H;0  molecules  held  by  these  strong  nonhydro- 
gen-bonded  forces  peak  just  above  3400  cm  -  V  Such  ex¬ 
tremely  low  NHB  frequencies  probably  occur  because  the 
O-H  bond  distance  lengthens,  and  the  stretching  force  con- 


gives  rise  to  a  frequency  just  above  the  centroid  of  the  NHB 
OH-stretching  component,  which  is  3628  cm  for  HDO. 
The  temperature  dependence  of  the  x-ray  data  is  also  consis¬ 
tent  with  the  observation  that  the  relative  NHB  OH-stretch- 
mg  intensity  from  HDO  increases  with  temperature  rise,67 
i.e..  that  the  concentrations  of  OH  oscillators  involved  in 
long  0-0  distnaces,  -  3.2  A.  increases  with  increasing  tem- 
perture. 

IV.  SUMMARY 

BF  correction  of  low-frequency  Raman  spectra  from 
water  and  thermodynamic  analysis  of  the  resulting  integrat¬ 
ed  intensities  yields  a  A//”  value  of  2.6  ±  0. 1  kcal/mol  OHO 
for  the  rupture  of  hydrogen  bonds.  This  AH  °  leads  to  reason¬ 
able  values  for  the  heat  capacity  and  for  various  entropies, 
and  is  in  agreement  with  XH°  values  obtained  from  vibra¬ 
tional  spectroscopy  for  the  one-  and  two-phonon  OH- 
stretching  regions,  as  well  as  from  other  methods  such  as 
depolarized  Rayleigh  scattering  or  molecular  dynamics. 
Some  implications  of  the  new  result  are:  that  the  low-fre¬ 
quency  intermolecular  Raman  spectrum  arises  predomi¬ 
nantly  from  partially  covalent  linear  and/or  weakly  bent  O- 
H  O  hydrogen  bonds.  HB;  that  the  breakdown  of  such 
partially  covalent  interactions  is  responsible  for  the  high 
heat  capacity,  compared  to  ice  or  to  steam;  and  that  this 
breakdown  produces  another  class  of  cohesive  interactions, 
NHB.  whose  further  complete  separation  would  require  an 
additional  3.2  kcal/mol  0-0  pair. 

Quantitative  (absolute)  Raman  spectra  obtained  in  the 
OH-stretching  region  yield  a  precise  isosbestic  point  at  3425 
cm"1.  Such  a  point  is  indicative  of  an  equilibrium  between 
the  HB  and  NHB  classes  of  interaction.  Raman  difference 
spectra  in  the  OH-stretching  region  show  a  rise  in  the  total 
integrated  intensity  above  the  isosbestic  frequency,  and  a 
decrease  below.  A  plot  of  the  coresponding  integrals  indi¬ 
cates  direct  proportionality.  This  fact  establishes  the  pres¬ 
ence  of  the  HB— NHB  equilibrium,  and  indicates  that  the 
OH-oscillator  concentration  is  conserved  in  terms  of  the  HB 
and  NHB  classes. 

V.  ADDENDUM 

A.  Volume  changes  and  pressure  effects 

The  equation 
In  [  /  /  f  28  1 5  —  /)] 

=  In  [  /«/ ( 1  —  fB  )  ] 

=  13.1 1 1  024 92"  -  3.126419  25, 

see  Fig.  2.  allows  values  of  7or/a  to  be  calculated.  Values  of 
fa  are  useful  in  making  rough  estimates  estimates  of  the  vol¬ 
ume  change,  AF  =  VB  —  F„,  corresponding  to  the  transfor¬ 
mation  U— B. 

The  volume  change  AFmay  be  approximated  from  the 
ideal  solution  formulation 

K  =/b  Yb  +/.  K  =faVB  +  (  I  ~fa )  (2D 

which  m  terms  of  the  volume  change  A  V  is 

V,  =f„dV+  F„.  (22) 


TABLE  III  Values  of  the  partially  covalent  hydrogen-bonded.  O-H  ■■  O. 
mole  fraction,  /* ,  at  various  temperatures  for  water. 


i.'C 

/. 

0 

0.842 

to 

0.818 

20 

0  794 

25 

0.781 

30 

0.768 

40 

0.743 

50 

0.717 

60 

0,692 

70 

0.667 

80 

0643 

90 

0619 

too 

0.596 

If  values  of fB  and  F;  are  known  for  widely  spaced  tempera¬ 
tures.  simultaneous  solutions  of  Eqs.  (21 )  or  (22)  yield  AF 
For  temperatures  of  4  and  90  °C.  which  span  the  present  low- 
frequency  Raman  data,  simultaneous  solution  yields  an 
average  AF  value  of  —  1.5  cmVmol  OHO.  Alternatively, 
Eq.  (22)  may  be  differentiated  with  respect  to  T  or  to  P 
Then  from  the  forced  assumption  that  all  volumes  except  F, 
are  independent  of  T  or  of  P,  and  from  substitution  of  the 
appropriate  expression  for  (dfB/dT)Por  (dfB  /dP)  T,  see  the 
derivation  in  Sec.  Ill  B  3.  values  of  A  F  vs.  t  may  be  calculat¬ 
ed  from  known  values  of  (dVi/dT) p  or  (dV,/dP)r.  When 
(dV,/dT)P  is  equated  to  A  FA Hfgfu/RT'-,  the  values  of  AF 
that  result  are  as  follows:  (1)0  *C,  +  0.26  cm'/mol  OHO; 
( 2 )  20  °C,  -  0.75:  (3)  40  *C,  -  1.4;  (4)  60  °C.  -1.9;  and 
(5)  80  °C,  —  2.5.  Of  course,  because  these  A  F  values  change 
both  in  sign  and  in  magnitude  from  0  to  80  °C.  it  is  evident 
that  the  assumption  relative  to  temperature  independence  of 
AF,  VB,  or  Vu  is  faulty,  but  all  of  the  calculated  A  F  values 
are  at  least  small.  Moreover,  the  CPC  result  agrees  with  a 
value  reported  by  Angel!.  Ref.  3b.  The  important  point  here 
is  that  if  AF  is  small,  its  magnitude  is  in  qualitative  agree¬ 
ment  with  Ref.  80.  where  no  significant  change  in  the  NHB 
OD-stretchtng  component  intensity  relative  to  the  HB  OD- 
stretching  component  intensity  could  be  ascertained,  within 
the  accuracy  of  the  computer  analysis  employed,  to  roughly 
10  kbar  at  32  °C. 

For  convenience  values  fo  fB  in  10°  increments  between 
0  and  100  °C  are  listed  in  Table  III. 

It  has  recently  been  possible  to  obtain  Raman  spectra  in 
the  fundamental  OH-stretching  region  from  strongly 
shocked  water  at  pressures  of  257  kbar  and  temperatures  to 
1437  °C  ( 1710  1C),  where  the  density  is  2  g  cm"’.'5  HB  and 
NHB  classes  or  interactions  are  both  definitely  indicated  up 
to  1437  °C  in  agreement  with  the  present  work.  However, 
most  of  the  hydrogen  bonds  are  broken  at  1437  *C,  and  the 
H-O  molecules  are  thought  to  be  held  together  by  other 
strong  forces,  e.g.,  electrostatic  (multipole),  induction,  and 
dispersion.  At  a  density  of  2  g  cm~\  the  OH-stretching  vi¬ 
brations  of  H-O  molecules  held  by  these  strong  nonhydro- 
gen-bonded  forces  peak  just  above  3400  cm"1  Such  ex¬ 
tremely  low  NHB  frequencies  probably  occur  because  the 
O-H  bond  distance  lengthens,  and  the  stretching  force  con- 


slant  drops,  as  the  0-0  nearest-neighbor  distances  in  the 
severely  bent  O-H  O  configurations  decrease.81 

B.  Polarization  measurements 

The  170  cm-1  Raman  band  from  water  has  been  as¬ 
signed  to  restricted  translation  of  H;0  molecules  along  the 
O-H  O  direction. The  60  cm-1  Raman  band  may  arise 
from  bending  of  the  O-H  —  O  units.  For  the  case  of  linear 
O-H  -  O  interactions,  point  group  Cx  the  H  •••  O  stretch¬ 
ing  vibration  v,X * ,  should  be  polarized,  p,  <  3/4.  v:2- 
should  be  depolarized,  p,  =  3/4.  However,  the  H  —  O  bond 
is  only  weakly  covalent,  and  thus  only  very  weak  polariza¬ 
tion  would  be  expected  for  v,2  +  .  This  weak  polarization 
may  be  difficult  to  detect  experimentally. 

Measurements  of  the  depolarization  spectrum,  i.e..p,  vs 
Av  were  made  in  this  work.  Broad  minima  were  observed 
near  Av  =  125  cm-'  and  Av  =  425  cm-1  (very  broad).  At 
125  cm-1,  p,  =0.71  ±  0.03  (very  weakly  polarized),  and  at 
425  cm-1,  p,  =  0.61  ±  0.04  (weakly  polarized).  The  fre¬ 
quency  of  425  cm  - 1  corresponds  to  the  libration  of  the  H-O 
molecule  around  its  two-fold  axis.17  At  the  O-H  —  O  bend¬ 
ing  frequency,  50  cm-1  (see  Fig.  3)  p,  =  0.74  +  0.01  (de¬ 
polarized).  Also,  a  sharp  minimum  in  p ,  was  seen  at  1650 
cm-1  which  corresponds  to  vyf,  (polarized).  Tor  p,  values 
in  the  OH-stretching  region  see  Ref.  13  and  especially  Ref. 
82. 
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